Spatial separation of state- and size-selected neutral clusters 
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We demonstrate the spatial separation of the prototypical indole(H20) clusters from the various 
species present in the supersonic expansion of mixtures of indole and water. The major molecular 
constituents of the resulting molecular beam are H2O, indole, indole(H20), and indole(H20)2- It is 
a priori not clear whether such floppy systems are amenable to strong manipulation using electric 
fields. Here, we have exploited the cold supersonic molecular beam and the electrostatic deflector to 
separate indole (H2O) from the other molecular species as well as the helium seed gas. The experi- 
mental results are quantitatively explained by trajectory simulations, which also demonstrate that 
the quantum-state selectivity of the process leads to samples of indole (H2O) in low- lying rotational 
states. The prepared clean samples of indole(H2 0) are ideally suited for investigations of the stereo- 
dynamics of this complex system, including time-resolved half-collision and diffraction experiments 
of fixed-in-space clusters. Our findings clearly demonstrate that the hydrogen bonded indole (H2O) 
complex behaves as a rigid molecule under our experimental conditions and that it can be strongly 
deflected. 

PACS numbers: 37.20. +j, 36.40.-c, 33.15.-e, 06.60.Ei 



The investigation of atomic and molecular aggregates 
- so called clusters - provides detailed information for 
bridging the knowledge gap between microscopic atomic 
and molecular systems and macroscopic samples. This 
includes, for example, the observation of ferroelectric- 
ity and superconductivity [1], the investigation of proton 
wires [2 J , and the evidence for single water molecules act- 
ing as catalysts of chemical reactions [3J . Acid dissolution 
has been studied at the microscopic level by stepwise at- 
tachment of solvents to molecules one by one [4J, and 
hydration shells are found to serve as heat sinks and to 
enhance the photostability of DNA |5J. Water clusters 
can have a significant effect on the chemistry and the 
climate in earth's atmosphere |6J. 

The actual dynamics of these processes could be inves- 
tigated, for instance, using ultrashort-pulse electron [7J|8] 
or X-ray diffraction [9], diffraction- from- within [10], or 
high-harmonic-generation experiments [llj. However, 
these experiments rely on clean, homogeneous samples in 
order to produce so called "molecular movies" for individ- 
ual cluster sizes and isomers. Therefore, it is highly desir- 
able to create pure samples of individual cluster species. 

For charged particles, size-selection can be performed 
according to the mass-to-charge ratio of the clusters. 
Large ionic systems can also be separated according to 
their shape in ion-mobility spectrometers [12J. The full 
separation of structural isomers (conformers) of specific 
sizes has so far been elusive, though. The preparation of 
state- and conformer-selected samples of (small) neutral 
molecules has a long history going back to the seminal 
contributions by Stern and Gerlach and by Rabi [T3HT5] , 
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For small molecules, static electric fields have extensively 
been used to select individual quantum states [T6] . More 
recently, switched electric and magnetic fields have been 
employed in the selection and slowing of individual or 
small sets of eigenstates of small molecules [17-20J. For 
large molecules, dynamic (strong) focusing has to be em- 
ployed pi| 122] . The rotational state-selectivity of this 
technique has been demonstrated [23 , 24J and it has been 
used for the spatial separation of individual conform- 
ers [25]. The state [26] and conformer [27] selection of 
large molecules has also been demonstrated using static 
inhomogeneous electric fields, going back to an original 
proposal for state selection by Stern [28J. 

It is a priori not clear, whether these techniques can 
be extended to complex molecules and clusters with large 
numbers of internal degrees of freedom, many of which 
have low excitation energies and will be excited even 
at low temperatures. Due to their large size (and cor- 
respondingly large moments of inertia) and their Sop- 
piness these systems have very high densities of states 
and, correspondingly, complex Stark energy curves with 
many real and avoided crossings. Indeed, it was con- 
cluded before that the "chaotic rotation" of such complex 
asymmetric rotors makes their manipulation very diffi- 
cult [29]. Here, we set out to demonstrate that the elec- 
tric deflection method is capable of spatially separating 
neutral polar molecules and clusters according to their 
size and conformation, providing unprecedented possi- 
bilities to investigate the structures and dynamics, even 
of complex aggregate systems. 

We have chosen to demonstrate the spatial separation 
of the prototypical indole(H 2 0) cluster from the vari- 
ous species in a supersonic cluster expansion of indole 
(CgHyN) and water (H2O) seeded in helium. Indole is 
the chromophore of the essential amino acid tryptophan 
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FIG. 1. (Color online): The experimental setup consists of 
a pulsed molecular beam source, a deflector, a time-of-night 
mass spectrometer and a pulsed dye laser. The inset shows 
a cut through the deflector and a contour-plot of the electric 
field strength for an applied voltage of 26 kV. 



and its intrinsic properties have been widely studied. Its 
emmission properties are regularly used in fluorescence 
studies of proteins, where spectral shifts are directly re- 
lated to the chromophores' environment [30 J. The influ- 
ence of water solvation on the indole chromophore has 
been discussed extensively, since it has strong influence 
on the electronic states, including the lowest-energy elec- 
tronically excited states, which are considered to inter- 
change their energetic order with the addition of water 
and other polar molecules |31| . The structure of the 
indole (H2O) complex was determined in the ground and 
electronically excited states [32j [33] and the influence of 
water on the lowest electronic states discussed [34H36] . 
Indole(H20)2 clusters were also investigated using so- 
phisticated multiple-resonance spectroscopic techniques 
that showed, for instance, water-to-7r-orbital binding mo- 
tifs [37J. The absorption spectra of indole in water solu- 
tions were interpreted using strong shifts and reordering 
of the molecules' electronic states |31j . 

Here, we experimentally demonstrate the preparation 
of a pure sample of the indole(H20) clusters using the 
electric deflector to separate this species from the clus- 
ter expansion soup. We show simulations that repro- 
duce the obtained results and will allow predictions for 
the preparation of further single-species samples even for 
such complex quantum systems. 

The experimental setup is shown in |Fig. 1| A super- 
sonic expansion of a few mbar of indole and a trace of 
water molecules seeded in 50 bar of helium is provided 
by a pulsed Even-Lavie valve [38]. The valve is operated 
at a temperature of 40 °C and at a repetition rate of 
20 Hz. Due to three-body collisions in the early phase of 
the expansion indole(H2 0) n clusters are formed. Under 
the employed conditions the indole(H20)2 signal (den- 
sity) is about an order of magnitude smaller than that 
of indole(H 2 0). Larger clusters are virtually absent from 
the beam. The resulting cluster beam passes a skim- 
mer (Beam Dynamics model 50.8) with an orifice diam- 
eter of 2 mm placed 22 cm downstream from the noz- 
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FIG. 2. (Color online): Structures, molecular constants, and 
Stark energies of (a) H2O, (b) indole, (c) indole(H20), and 
(d) indole(H2 0)2; see text for details. 



zle and enters a differentially pumped vacuum chamber 
which houses the electric deflector. The length of the 
deflection device is 24 cm and the vertical distance be- 
tween the two electrodes, a stainless steel trough and 
rod, is 2.3 mm. We typically apply voltage differences 
of kV to 26 kV between the two electrodes. This re- 
sults in electric field gradients as shown in Fig. 1, with 
field strengths and gradients of up to 120 kV/cm and 
250 kV/cm 2 in the center of the selector, respectively. 



The electric field gradient exerts a force on the molecules 
and clusters by which they are deflected towards the rod 
if they are in high- field-seeking states. Since the ori- 
fice diameter of the first skimmer and the distance of 
the two electrodes are comparable, we have moved this 
skimmer slightly towards the trough. This avoids the 
deflected beam crashing into the rod electrode and im- 
proves the sensitivity to deflected molecules. After pass- 
ing a second skimmer with an orifice diameter of 2 mm 
mounted on a motorized translation stage the molecules 
enter the detection chamber. This chamber houses a 
time-of-flight (TOF) mass spectrometer. The position of 
the second skimmer is adjusted to maximize the trans- 
mission through the skimmer when the molecules are de- 
flected. The molecular beam is crossed by a laser at 
90°, in the center of the TOF-mass spectrometer. In- 
dole and indole-water clusters are ionized via one-color 
resonance-enhanced multi-photon ionization (REMPI) at 
around 283 nm exploiting the lowest allowed electronic 
excitation of the indole chromophore. Due to the reso- 
nant excitation step involving different laser frequencies, 
the ionization is spectrally fully species selective. The 
laser light is provided by a nanosecond Nd:YAG laser 
(Innolas) pumped frequency doubled dye laser (Fine Ad- 
justment). The light is focussed by a 75 cm-focal- length 
lens and the typical pulse energy is a few mJ. For detec- 
tion, the resulting ions are accelerated towards a multi- 
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channel plate (MCP) detector and time-of-flight (TOF) 
mass spectra are recorded. Indole densities are deter- 
mined from mass 117 (indole + ) and indole- water den- 
sities are determined by summing masses 117 and 135 
(indole + and indole(H 2 0) + ) following (neutral- species- 
selective) resonance-enhanced ionization of the respec- 
tive molecules or clusters. The focusing lens is mounted 
on a translation stage and spatial density profiles of the 
molecular beam are obtained by moving the lens perpen- 
dicular to the flight direction (^-coordinate, see Fig. 1) 
of the molecular beam. 



Fig. 2| shows structures, molecular constants, and Stark 
energies of all relevant species. These are experimen- 
tally determined for H 2 [39l HO], indole [36l SU, and 
indole(H 2 0) [33, 36J and based on ab initio calculations 
(GAMESS-US [42], B3LYP/6-31+G*) for indole(H 2 0) 2 ; 
the calculated structure and dipole moment are the same 
as the ones reported by Carney et al. [43, 44J. The avail- 
able experimental quartic centrifugal distortion constants 
for water, indole, and indole (H 2 0) are used in the cal- 
culations [45 ] l46 j . However, the resulting calculated de- 
flection profiles using the rigid-rotor (not shown) and the 
semi-rigid-rotor approximations are identical. The inset 



of Fig. 2 a is a zoom into the Stark energies of water in 
order to visualize its small Stark shifts. All populated 
states of the molecules and clusters, at our experimen- 
tal conditions, show a decrease of energy with increasing 
field strength for the relevant range of ~50-150 kV/cm. 
Thus, all molecules are in high-field seeking states. The 
change in energy in the strong electric field depends on 
the specific rotational quantum state. In general, the 
Stark effect decreases with increasing J quantum number 
giving rise to the largest change in energy for the rota- 
tional ground states. Therefore, molecules in the lowest 
rotational states are deflected the most. Indole (H 2 0) 
shows the largest Stark effect of all molecules and clus- 
ters under consideration - its potential energy changes 
50 % more than for indole (H 2 0) 2 and twice as much as 
for indole: If the clusters behave "well" (vide supra), the 
biggest change in the molecular beam profile, due to the 
interaction with the inhomogeneous electric field in the 
detector, clearly is expected for indole(H 2 0). 



Fig. 3| shows the indole and indole(H 2 0) density pro- 



files of the molecular beam, obtained by scanning the fo- 
cused REMPI laser beam perpendicular to the molecular 
beam direction, for deflector voltages of 0, 10 kV, 20 kV 
and 26 kV. The width (FWHM) of the undeflected beam 
is about 1 mm, limited by the mechanical aperture be- 
tween the first skimmer and the deflector trough in order 
to avoid collisions of the molecular beam with the rod and 
to increase the sensitivity to detect deflected molecules. 
The density profiles show that all molecules and clusters 
are deflected upwards and that they are dispersed as the 
deflector voltage is increased. This can be attributed to 
all quantum states being high-field seeking for electric 
field strengths present in the deflector and a variation 
(decrease) of the Stark effects with increasing rotational 
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FIG. 3. (Color online): Deflection curves of (a) indole and 
(b) indole(H20) for deflector voltages of (green diamonds), 
10 kV (yellow circles), 20 kV (blue triangles) and 26 kV (red 
squares). The lines show the outcome of the respective tra- 
jectory calculations; see text for details. 
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FIG. 4. (Color online): The torsion potential and energy level 
scheme of indole(H20) in its vibronic ground state (excluding 
torsion) . 



respond to simulated molecular beam profiles [46J. The 
force field inside the deflector has been derived from the 
calculated electric field and the Stark curves and the pro- 
files have been obtained via Monte-Carlo trajectory sim- 
ulations of individual molecules. The good agreement 
between simulation and experimental data has been ob- 
tained by assuming a rotational temperature of 4 K for 
indole and of 6 K for indole (H 2 0) and by positioning the 
source and the two skimmers accordant to the experi- 
mental conditions. 

The indole and water moieties of indole(H 2 0) are con- 
nected by a single hydrogen bond, and they can rotate 
with respect to one another about this single bond. This 
torsional motion introduces the complication that the 
cluster generally cannot be treated as a (semi)rigid ro- 
tor [47J. However, under our experimental conditions 
the cluster regains rigidity due to the low temperature 
of our sample, the large rotational constant of the in- 
ternal rotation of water, and the (dis) allowed symmetry 
couplings of the populated states: Consider the twofold 
symmetric potential energy curve of indole (H 2 0) shown 
The angle between the molecular plane of in- 
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excitation, respectively; see Fig. 2 The lines in Fig. 3 cor- 



dole and the molecular plane of the water moiety is de- 
noted by a. The torsional energy levels are labeled by 
the torsional quantum number v t . Each torsional level 
is split due to tunneling splitting and the corresponding 
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torsional sublevels are denoted by a = 0,1. The barrier 
height, ^2=168. 5 cm -1 , the angle of the axis of rota- 
tion of the water moiety in the be symmetry plane of 
the cluster (corresponding to the ab plane of bare indole) 
with respect to the b axis, /3=55°, as well as the rota- 
tional constants have been taken from reference [32j In 
order to determine the energy levels we have numerically 
solved the one-dimensional Schrodinger equation [481149], 
The included assumption that the rotating moiety is a 
symmetric top molecule is not strictly fulfilled for water. 
However, the principal axes of indole(H20) do not change 
significantly when the water rotates and we can treat wa- 
ter as a symmetric top [32j[50]. With this model, we cal- 
culate a subtorsional splitting in v t = of 0.30 cm -1 in 
agreement with the originally derived value of 0.31 cm -1 
for this model [32J. In order to estimate the contribution 
of v t > 1 due to its initial state population prior to en- 
tering the strong field, we assume the population of the 
energy levels to follow a Boltzmann distribution. The en- 
ergy difference between \v t ^cr) = |0,0) and \vt,<j) = |1, 1) 
is 70.9 cm -1 and the difference between \v t: cr) = |0, 1) 
and \vt,&) = 1 1,0) is 76.5 cm -1 . This gives rise to a rel- 
ative population of v t = 1 with respect to v t = of less 
than 1 % at vibrational temperatures at 20 K or below, 
as is expected for our supersonic expansion. Therefore, 
we neglect population with v t > 1 in our calculation, 
confirmed by the good agreement with the experimental 
results. 

The two subtorsional states in v t = do not couple 
in the strong electric field of the deflector due to their 
symmetry (see appendix [A|. Coupling to v t > 1 can be 
neglected, because the field in the deflector is not strong 
enough to efficiently couple v t = to these states (es- 
sentially because the energetic separation is too large). 
Therefore, under our well-controlled experimental con- 
ditions indole(H2 0) can be treated as a (quasi) rigid 
molecule. 
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~5| shows the experimental and simulated 
beam deflection profiles for indole, indole(H20), and 
indole(H20)2, and the simulated profiles for water, for 
voltages of and of 26 kV applied across the deflector. 
The deflection-field-free density profiles are depicted in 
the inset of|Fig. 5[ All undeflected and deflected beam 



profiles have been normalized consistently to the field- 
free peak intensity for each species. This accounts for 
the unknown (relative) cross sections for ionizing the 
molecules and clusters and for the unknown number den- 
sities of the different species. The width of the unde- 
flected beams is the same for all species. For strong 
electric fields applied, the beam profiles of all species 
broaden - due to the different effective dipole moments 
of the populated states - and shift upwards to larger y 
values. Indole (H2O) is clearly more deflected than wa- 
ter, indole, and indole(H20)2, as it is expected from the 
Stark curves shown in |Fig. 2[ We note that "exception- 
ally strong deflection" was reported for water in two of 
its J = 1 states [51J. The reported behavior is repro- 
duced by our simulations, but this level of deflection is 
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FIG. 5. (Color online): Normalized experimental (symbols) 
and simulated (lines) molecular beam profiles for water (dot- 
ted line), indole (blue circles/dashed line), indole(H20)(red 
squares/solid line), and indole (H2 0)2 (green triangles/dash- 
dotted line) at a deflector voltage of 26 kV. The inset shows 
the corresponding undeflected normalized beam profiles. 



not very strong on the scale of our experiment. These re- 
sults, i.e., the profiles shown in|Fig. 5| demonstrate that 



it is possible to perform experiments on pure samples of 
indole (H2O) by simply selecting the species present in the 
grange from 2.3-3.3 mm. At y = 2.75 mm (3.0 mm) the 
signal is 10 % (4 %) of the undeflected peak intensity with 
only the 290 (110) lowest rotational quantum states pop- 
ulated, instead of the 4600 states of indole(H20) popu- 
lated in the original molecular beam. Such state-selected 
samples have been shown to be very amenable to align- 
ment and orientation control [26| 152]. 

In conclusion, we have spatially separated indole (H2O) 
from a cluster beam containing water, indole, and var- 
ious indole- water clusters as well as atomic helium. 
We have experimentally demonstrated that indole (H2O) 
is deflected more than the indole monomer and the 
indole(H 2 0)2 cluster. We have quantitatively simulated 
all experimental data using trajectory simulations. The 
good agreement of the simulations with the experiments 
demonstrate that the manipulation of complex molecu- 
lar systems and even weakly bound clusters is possible 
and can be understood - if the initial temperature of the 
sample is low enough. These simulations also demon- 
strate that the indole(H2 0) clusters are separated from 
H 2 and, obviously, from the (undeflected) atomic he- 
lium seed gas as well as from larger water, indole, or 
indole- water clusters, whose contribution to the beam is 
minimal and which have smaller dipole-moment-to-mass 
ratios. Due to the state-specific deflection, it will also be 
possible to separate different structural isomers, which 
will occur for larger clusters [43] - analogous to the con- 
former separation of individual molecules [25j [27] . 

Recently, in a related experiment the selective 
alternating-gradient focusing of benzonitrile and the 
benzonitrile-argon cluster was demonstrated [53J. 
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The resulting pure samples of individual cluster species 
will allow for novel experiments investigating the struc- 
tures and dynamics of these species. One particular ex- 
periment we foresee is the investigation of the energet- 
ics and stereodynamics, and their correlation, of half- 
collisions of indole(H20) in the molecular frame. One 
would 3D-orient [52] the molecules in the deflected pure 
sample, photodissociate the cluster species using various 
photon energies, and, successively, ionize the neutral wa- 
ter fragment from the reaction. The full momentum vec- 
tor of the water molecule in the molecular frame would 
then be determined using ion imaging. Similar investi- 
gations using time-resolved gas-phase X-ray and electron 
diffraction experiments are also envisioned [54]. Such 
experiments would provide novel details about the in- 
termolecular interactions during, and their influence on, 
chemical reactions of complex molecules [55H57] . 
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Appendix A: Symmetry properties of indole(H2 0) 
and coupling between sub-torsional levels 

Here, we consider the symmetry properties of the tor- 
sional and sub-torsional levels of indole (H2O) and the 
couplings between them due to an external dc electric 
field. The argumentation is equivalent to the case of 
phenol(H20) [58]. Because it is important regarding the 
Soppiness or rigidity of the indole(H2 0)-cluster in our 
experiment the derivation is provided here in detail. 

We describe the internal motion of the water moiety 
with respect to the indole moiety using the principal axis 
method (PAM) [48j [49| page 582]: We use the principal 
axes of the cluster as a coordinate system in combination 
with an internally rotating symmetric top moiety. The 
axis of internal rotation coincides with the symmetry axis 
of that symmetric top and, therefore, internal rotation 
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TABLE II. The symmetry classification of eigenstates of field- 
free indole(H 2 0) and the respective nuclear spin statistical 
weights. 



does not change the moments of inertia of the overall ro- 
tation nor the cluster's principal axes. This condition is 
strictly speaking not fulfilled for the case of indole (H2O) 
since the rotating water moiety is an asymmetric top 
molecule. However, the moment of inertia of the rotating 
water molecule about its pseudo-symmetry axis is small 
compared to the moment of inertia of the total molecule. 
The principal axes and moments of inertia of the clus- 
ter are only slightly and insignificantly altered by the 
rotation of the water molecule. The applicability of this 
approach is confirmed by related spectroscopic studies of 
hydrogen bound molecular clusters [32 , 36l l50| 158 ] . 

In field-free space, indole(H20) belongs to the molecu- 
lar symmetry (MS) group G 4 = {E, (12), E* , (12)*} [59J 
and the corresponding character table is provided in |Ta-| 
ble I The symmetry of the total wave function T tot is of 



B\ or B2 symmetry since the total wave function must 
be antisymmetric with respect to the exchange of the 
two protons (operation (12)). The symmetry of the elec- 
tronic, vibrational, torsional, rotational and nuclear spin 
states are summarized in |Table II] In the last column, 



the resulting nuclear spin statistical weights (nssw) are 
given. 

In an electric field the molecular symmetry group of 
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TABLE I. Character tables for the molecular symmetry 
groups G 4 and G 2 [58] [59] 



TABLE III. The symmetry classification of eigenstates of 
indole(H 2 0) in a dc electric field and the respective nuclear 
spin statistical weights. 
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indole(H 2 0) reduces to G2 = {£",(12)}, which contains 
only two irreducible representations, A and B; the char- 
acter table is given in ITable II 
G2, 



The symmetries, under 



of all rotational quantum states of indole (H2O) in 
its vibronic ground state are obtained directly by corre- 
lating Al } 2 and B12 symmetries in G4 to A and B in 
G2, respectively. They are provided in |Table III| Tor- 
sional sublevels with the same quantum number a have 
the same symmetry (A for all a = sublevels and B for 



all a = 1 levels). The nonzero \i a and /i& components of 
the electric dipole moment are both of A symmetry, be- 
cause they do not change sign under the operation (12), 
i. e., (12) does not change the dipole moment of the clus- 
ter. Thus, an electric field mixes only states having the 
same a. Consequently the sub torsional states within a 
given v t do not mix. In our experiment more than 99 % 
of the population resides in the v t = states, see |page^} 
and indole(H 2 0) can be treated as a rigid rotor molecule 
under our experimental conditions. 
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